Inhomogeneous spacetimes as a dark energy model 



David Garfinkle 

Dept. of Physics, Oakland University, Rochester, MI 48309, USA 
E-mail: garf inklOoakland . edu 

Abstract. Lemaitre-Tolman-Bondi inhomogeneous spacetimes are used as a 
cosmological model for type la supernova data. It is found that with certain parameter 
choices the model fits the data as well as the standard ACDM cosmology does. 
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1. Introduction 

Observations of type la supernovaejH Ej have led to the conclusion that the expansion 
of the universe is accelerating. Since any matter with non-negative pressure causes a 
deceleration, this leads to the conclusion that the universe contains a substantial amount 
of an exotic type of matter refered to as "dark energy." A cosmological constant can 
account for the requisite dark energy, and indeed a model with a cosmological constant 
as well as ordinary baryonic matter and cold dark matter, the ACDM model, has become 
the standard cosmology. Though it agrees quite well with the data, the ACDM model 
has certain disturbing features: From the point of view of Planck units, the cosmological 
constant seems unnaturally small by about 120 orders of magnitude. In addition, the 
ACDM model requires that we live in the very special cosmological era when matter and 
dark energy have comparable densities. These disturbing features motivate the search 
for alternatives to the standard ACDM model. One such alternative is the theory 
of Kolb et a/[Sl EJ E] that the universe contains only ordinary matter and that the 
apparent acceleration is an effect of the universe's inhomogeneities. However, the work 
of j3 IH Ej has been criticized on many grounds. [HI 13 Ej In particular, their treatment 
uses perturbation theory beyond its regime of validity since the effects that they try 
to account for are large. Since a full non-perturbative treatment of an inhomogeneous, 
anisotropic universe is quite difficult, some authors [01 EHJ HU EE21 EH IHj have treated 
a toy problem: an inhomogeneous but spherically symmetric universe of pressureless 
matter. These models, the Lemaitre-Tolman-Bondi (LTB) spacetimes. |15 [ ITT) | IT7j are 
well studied, so all that is needed is to apply them to the question of the apparent 
acceleration of the universe. This has been done by several authors: P ITUl ITU H2^ 1131 IT4] 
however, for the most part the authors of these references have used purely formal 
notions of the acceleration of the universe. In contrast, Vanderveld et al lSj and Alnes 
et al[19\ have emphasized that it is more appropriate to see whether the model fits the 
relation between redshift and luminosity from which the cosmological acceleration was 
deduced. 

In this paper, we will treat a class of LTB spacetimes and calculate their 
compatiblity with the observational data ofpQ. Section 2 contains a presentation of 
the notation and methods to be used. Results are presented in section 3 and their 
implications discussed in section 4. 

2. Luminosity distance in LTB models 

The LTB metric takes the form [20 

ds 2 = -dt 2 + Y^rjdr 2 + r 2 (d6 2 + sm 2 6d<j) 2 ) (1) 

Here the area radius r is a function of both coordinates t and f and a prime denotes 
derivative with respect to f. The function / depends only on f . It follows from Einstein's 
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equation that 

r 2 = f + - (2) 
r 

where an overdot denotes derivative with respect to t and where the function F depends 
only on f. The density is given by 
F' 

P = (3) 

It is helpful to introduce the quantities a, A and B by 

r = ar (4) 

/ = Ar 2 (5) 

F = Bf 3 (6) 

Then equation (J2J) becomes 

a 2 = A + — (7) 
a 

The solution of equation (|ZJ) is 

f a du 



where to is a function of f whose physical meaning is the time at which the shell of dust 
labeled by f has zero radius. 

At first sight, it seems that the LTB solutions depend on three arbitrary functions: 
A, B and to- However, there is still the coordinate freedom to change f to any function 
of f. We will use that coordinate freedom to make B constant. This makes the 
solutions with constant A and zero t Q the standard Friedman- Robert son- Walker (FRW) 
homogeneous and isotropic cosmologies. We will further choose the value of B to be 4/9 
so that the A = FRW cosmology takes the standard form a = t 2 ! 3 . In addition we will 
consider only models with to — 0. These models have a genuine Big Bang singularity in 
the sense that all dust shells are at zero radius at the same time. 

Define the null vector l a by 

Then a future directed null geodesic takes the form k a = —uoi a and the redshift z is 
given by 1 + z = uj/uj where a zero subscript denotes the value of the quantity at f = 0. 
The luminosity distance is given by [TH] di = r(l + z) 2 '. 

We would like to find the luminosity distance as a function of redshift. To that 
end, it is helpful to find how each of these quantities varies as one goes along the past 
light cone of the observer. It follows from the geodesic equation that 

tV a z= r -{l + z) (10) 
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Since di is a function of z and r we need to know also how r varies along the past light 
cone. We have 



fV tt r=^l + /-r 
It is also helpful to know how a varies along the past light cone. We have 

tV a a = ^TJ^ - d 

Thus along the past light cone we find 
dr 

dz 
da 
dz 



[l + z)f' 
1 



1 + / 



l + fa' 



ar 



(11) 
(12) 

(13) 
(14) 



;i + z)f 

If we can find all the quantities on the right hand sides of equations (flUj) and (fT4^) 
as functions of z, r and a then we will be able to numerically integrate these equations 
and thus find 4 asa function of z. Since f = r/a we can find any function of f and thus 
given a choice for the function A we can find /. Similarly, we can find f using equation 
(J2J) and a using equation (J7J). Given a' and a' we have r' = a + fa' and f = a + fa' so 
all that remains then is to find a' and a'. Differentiating equation (jHJ) with respect to f 
and solving for a' we obtain 




Where / is given by 
, 1 



a A B 

2Jo V u 



-3/2 



du 



This integral can be done in closed form, and for A > takes on the value 



B 



2A 5 /2 



Aa 
~B 



In 



l + s 
1 - s 



where 



8=11 



B 

Aa 



-1/2 



(15) 



(16) 



;i7) 



(18) 



Now differentiating equation ([To]) with respect to t and using equations (0) and (fTB^l we 
obtain 



A' 
Y 



a 



-1/2 



(19) 



3. Results 

In [THj models with A = and to 7^ were considered, while jTH] treat models with both 
A and to nonvanishing. In contrast, we will consider models with to = and A ^ 0. 
There remains the choice of the function A. Observations of the Cosmic Microwave 
Background suggest that fl = 1 whereas observations of galaxy clusters suggest that Q 
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in matter is approximately 30%. In the context of the LTB models, this suggest that 
A should approach zero at large distances and be positive at smaller distances. We will 
choose A to have the form 

A = — (20) 
1 + (cr) 

where the constant c will be chosen for the best fit with the data. Let Qm be the value 
of Q in matter at our present time and position. Then in the LTB models with our 
choice of A and B we have 

fiM=(l + ^)" 1 (21) 
Thus, given a choice for the value of Qm we find that the current value of a is given by 

= ^(0^-1) (22) 

The initial values of r and z are zero. We choose the value of Qm and then use equation 
(J22J) to get the initial value of a. We then numerically integrate equations (fTSJl and (jUj) 
outwards to find the luminosity distance as a function of z. 

In P what is computed is the effective magnitude, which is 

m effective = Mg + 6 \ O g(H d L ) (23) 

Here H is the Hubble constant and M.b is a constant to be found by the best fit to 
an FRW model. In our model, Hq is simply the initial value of a/ a. In (T| the models 
are all FRW and the free parameters are Q\ and Qm- The best fit, the one with the 
smallest x 2 is f° r the model with Q^ = 1.32 and Qm = 0.73. For any other model, 
whether FRW or not, one can then use the supernova data of [T] to compute its \ 2 an d 
then subtract off the x 2 of that best FRW fit to obtain A% 2 which is a measure of the 
goodness of fit. In particular, the standard ACDM model has = 0.74 and Qm = 0.26 
and has A% 2 = 1.2. 

We will consider two different models: one with Qm — 0.3 and one with Qm = 0.2. 
In each case the value of c will be chosen for best fit with the supernova data. The errors 
in both m ^ ectlve and z are used to compute y 2 . In figure ((TJ) are plotted the supernova 
data and the Qm = 0.3 LTB model. For this model c = 8.5 and A% 2 = 3.8. Thus, the 
standard ACDM model (the one with Q^ — 0.74 and Qm = 0.26) is a better fit to the 
data. However, given the possiblity of systematic uncertainties in the supernova data, 
it not clear whether this difference is significant. In figure (J2J) the Qm = 0.3 LTB model 
and the standard ACDM model are plotted. Note that the curves agree very closely. 

In figure are plotted the supernova data and the Qm = 0.2 LTB model. For this 
model c = 5.1 and A% 2 = 0.3. This is an even better fit to the data than the standard 
ACDM model. In figure JU the Q M = 0.2 LTB model and the standard ACDM model 
are plotted. Note that the curves are almost identical. 

In jTH] it is emphasized that a physically reasonable LTB model must be smooth 
both at the origin and elsewhere. For the coordinates that we have chosen, smoothness 
at the origin is insured provided that A is a smooth function of f 2 . Smoothness elsewhere 
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Figure 2. Plot of effective magnitude vs redshift for the standard ACDM model (solid 
curve) and the = 0.3 LTB model (dashed curve). 
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Figure 4. Plot of effective magnitude vs redshift for the standard ACDM model (solid 
curve) and the Hm = 0.2 LTB model (dashed curve). 
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Figure 6. Plot of A vs rcdshift for the LTB models with = 0.3 (solid curve) and 
the Qm = 0.2 LTB model (dashed curve). 
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requires that there be no " shell-crossing" ie that r' not vanish. In figure (JHJ) are plotted 
r' vs redshift along the past light cone for the Qm — 0.3 and Qm — 0.2 LTB models. 
Since r' does not vanish, there is no shell crossing. 

The fairly large values of the constant c in these models means that the spatial 
size of the underdensity in which we live is fairly small. In figure (JBI) are plotted A vs 
redshift along the past light cone for the Qm — 0.3 and Qm — 0.2 LTB models. Note 
that A = corresponds to an Q = 1 FRW spacetime. Thus in these models, we inhabit 
an underdensity of fairly small spatial extent in what is otherwise a spatially flat FRW 
cosmology. 



4. Conclusions 



This work shows that an inhomogeneous but spherically symmetric cosmological model 
can account for the apparent acceleration in the supernova data without any dark energy. 
The models have a simple form for the inhomogeneity and they fit the supernova data 
as well as the standard FRW model with a cosmological constant does. 

However, this does not make the LTB spacetime a viable cosmological model. There 
is more cosmological data than the supernovae. In particular, data from the cosmic 
microwave background and from weak lensing indicate that the universe contains a large 
amount of matter that is "unclumped." This would be accounted for by a cosmological 
constant; but is unlikely to be accomodated by a LTB model or any simple modification 
of one. Furthermore, since in these simple models, we inhabit an underdensity of farily 
small spatial extent, the models may not be compatible with data for galaxy counts as 
a function of redshift. Nonetheless, the work ofJTH] suggests that it might be possible to 
accomodate some of this additional cosmological data within the LTB models. These 
issues needs further study. 
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